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Macroscopic description of microscopically strongly inhomogenous 

systems 

A mathematical basis for the synthesis of higher gradients metamaterials 


A. Carcaterra^’^’^, F. dell’Isola^’^, R. Esposito^ and M. Pulvirenti^’^ 


Abstract 

We consider the time evolution of a one dimensional n-gradient continuum. Our aim is to construct 
and analyze discrete approximations in terms of physically realizable mechanical systems, called mi¬ 
croscopic because they are living on a smaller space scale. We validate our construction by proving a 
convergence theorem of the microscopic system to the given continuum, as the scale parameter goes to 
zero. 


1 Introduction 

Continua with exotic behaviors are acquiring an increasing attention for their interest in technological 
applications (see e.g. [in [Ml [la E [Ml ES] and references therein). In this paper we address what, in a 
sense, is an inverse problem: given a continuum model we seek for those mechanical systems which, at a 
certain length scale, behave as specified by the chosen continuum model. The aim is to understand the 
microscopic properties of such systems to obtain information on how to realize (synthesize) them, at least 
in principle. 

To be more precise, we are interested in a metamaterial which, roughly speaking, is an array of 
elementary individuals, much smaller than the typical macroscopic size, arranged in periodic structures 
and exhibiting unusual macroscopic behavior. 

In our mathematical analysis we want to consider such a continuous system as described by a partial 
differential equation generated by a Lagrangian which summarizes all the macroscopic properties we may 
desire. Then we discretize this system and manage to identify such a discretization as a real conservative 
mechanical model. In other words we start from a macroscopic behavior and describe one possible 
microscopic interaction which realizes it at a macroscopic level. Finally we give a mathematical foundation 
to this procedure by proving a convergence result. 

From a mathematical point of view, we underline once more that this is an inverse problem, compared 
to the one (largely unsolved) formulated by D. Hilbert in his famous speach in 1900 at ICM in Paris (see 
m) in which he encouraged to prove rigorously the transition from particle systems to fluid dynamics 
(Hilbert’s 6-th problem). However it is worth to stress that we are working in the framework of continuum 
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mechanics, but our microscopic elements, even if small in macro unities, are large compared with molecular 
scales. 

We conclude this introduction by spending some more words on metamaterials, collocating them in 
the framework of generalized continua, with a particular emphasis to the pioneering work of G. Piola (see 

[911251 [ 2 ]). 

The rest of the paper is organized as follows. In Section 2 we introduce continuous and discrete 
Lagrangians and discuss the identihcation problem, namely we specify the mechanical systems outlined 
by the discretization procedure. In Section 3 we formulate and solve the associated convergence problem. 

We remark that our work concerns one-dimensional systems only. This is of course a severe limitation, 
but, on the other hand, it is a natural setting to start with. 

1.1 Mechanical metamaterials 

By suitably rephrasing Engheta and Ziolkowski m and Zouhdi et al. [29], metamaterials are materials 
which are first theoretically conceived and then engineered to have properties very unlikely to be found 
in nature. 

They are obtained by suitably assembling multiple individual elements constructed with already avail¬ 
able microscopic materials, but usually arranged in (quasi-)periodic sub-structures. Indeed the properties 
of metamaterials do not depend only on those of their component materials, but also on the topology of 
their connections and the nature of their mutual interaction forces. In literature it is currently specified 
a particular class of metamaterials, so called mechanical metamaterials, those in which the particular 
properties which are “designed” for the newly synthesized material are purely mechanical. The present 
paper deals exactly with such a class. 

We explicitly remark here that in the present paper we use the adjective “microscopic” or “micro-” 
meaning all those length scales which are (much) smaller than the scale at which continuum mechanics 
is applicable. In particular we do not attach any value in SI units to each considered length scale. 

The particular shape, geometry, size, orientation and arrangement of the elementary individual el¬ 
ements can affect, for instance, the propagation of waves of light or sound in a not-already-observed 
manner. In this way one can create material properties which cannot be found in conventional materials. 

Particularly promising are those micro-structures which present high-contrast in microscopic proper¬ 
ties. These structures, once homogenized, have shown to produce generalized continua (see e.g. umm)- 
These micro-structures, although remaining quasi-periodical, are conceived so that some of the physical 
properties which are characterizing their behavior are diverging when the size of the representative ele¬ 
mentary volume tends to zero, while simultaneously some others are vanishing in this limit. 

To give a hint of the possible applications of newly designed metamaterials we list here some among 
the papers which are more relevant to our results, especially in the perspective of their extension to 2D 
and 3D systems. In [TS] it is shown how to synthesize a composite medium exhibiting negative effective 
bulk modulus, negative effective mass density (see also 0), or both properties. In |16] materials with 
negative Poisson’s ratio (auxetics) was designed, and they were fabricated in 1999 (see Xu et al.[28] ). 
One of the most famous examples of such materials is the Goretex whose negative Poisson ratio opened 
unexpected possibility to e.g. vascular surgery. 

The damping effects can be also suitably designed using special selection of the material microstructure 
as reported in [am, or the acoustic and optical effects as negative refraction, lensing and cloaking mm- 

All described materials can be modeled at a micro-level as finite dimensional Lagrangian systems and 
their effective properties were all obtained via a kind of homogenization procedure. 


1.2 Generalized Continua 

In the first half of XIX century the design of structures became an intellectual activity based on the 
rigorous application of predictive mathematical models. These models were formulated by means of a 
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precise postulation process and originated a series of problems or exercises directly motivated by the 
engineering applications, which were solved by means of the use of the then newly developed techniques 
of mathematical analysis. 

The model describing the mechanical behaviour of materials introduced by Cauchy - although very 
accurate for a large class of phenomena - cannot be applied to all materials in every physical condition. 

More general models were formulated by Gabrio Piola in the same years, but only recently they were 
considered in engineering for applications. 

In some formulations of continuum mechanics, the possibility of the dependence of deformation energy 
on higher gradients of displacement, is rejected, due to an apparent (see 0) incompatibility with the 
second principle of thermodynamics ([10],[H]). On the other hand, physicists, for instance Landau m, 
always considered this dependence as admissible, as they are accustomed to base the postulation of 
physical theories on the principle of least action or on the principle of virtual works, which is exactly the 
same starting point of G. Piola [25]. 

Actually, when introducing Piola continua, the true conceptual frame settled by Cauchy, Navier and 
Poisson is to be drastically modified. The concept of stress becomes secondary and the main role is played 
by deformation measures together with action and dissipation functionals. The Euler-Lagrange equation 
obtained in this more encompassing modeling process cannot be anymore regarded to coincide with the 
balance of force unless one generalizes the concept of force. This can be done by introducing generalized 
actions as the dual quantities in the work of the gradients of displacements (see e.g. na Eoi [221 EH da 

EM)- 

Actually the same concept of contact interaction has to be completely modified, and the crucial point 
of determining the correct boundary conditions which can be assigned in generalized continua theory has 
been addressed only very recently (see e.g. following the original ideas by Piola [25]. 

2 Microscopic and Macroscopic descriptions 

In what follows we will consider two length scales I and L with I L. We will call microscopic or micro 
the description at the length scale I, while macroscopic or macro will be the attribute relative to the 
description which is suitable at the lenght scale L. 

We assume that the most suitable micro-description at micro-scale is “discrete” i.e. based on the model 
“material particle” (as done by Poisson, Navier and -in some works- by Piola), while the description which 
has to be used at the macro-level is that of a continuum, as introduced e.g. by Lagrange, Cauchy or again 
Piola. 

Remark however that we will not limit our attention to systems which verify the assumptions put 
forward by Cauchy and Navier. We will consider, actually, those continua which have been considered 
by Piola (and then by many others, including Toupin, Green, Rivlin and Mindlin) i.e. so called higher 
gradient continua. 

To quantify the above considerations we will introduce, in the sequel, a small parameter e > 0 
indicating the ratio between typical micro and macro scales, possibly to be sent to zero to outline a 
suitable asymptotic behavior. 

2.1 The basic macroscopic continuous model 

Let / = (0, L) C M be a finite interval assumed as reference configuration of the considered one-dimensional 
continuum. We label each element of the continuum with the coordinate a: € I of its placement in the 
reference configuration. The actual configuration of the continuum is described by the displacement field 
u = u(x, t) which represents the horizontal displacement at time t of the element x from its position in 
the reference configuration. 

Fixed an integer n > 1, for such a system we introduce the Lagrangian 
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C[u,u) = ^ J \u{x)\‘^ — J ^{u{x),Du{x),... ,D^u{x)). ( 2 . 1 ) 

Here, D^u is the /c-th x-derivative of u and 

( 2 . 2 ) 

is a function whose properties will be specified later on. 

Note that ^{u, Du,..., D'^u) is the potential energy density corresponding to the dispacement u and 
describes the constitutive properties of the medium under investigation. 

The action on the time interval (0, T) is consequently defined as 

A = [ C{u{- ,t)), (2.3) 

Jo 

where u{x, t) = dtu{x, t) is the time derivative. 

To deduce the Euler-Lagrange equations from the stationary action principle, we have first to specify 
the kinematic boundary condition for our problem. In the sequel we shall assume either 

• periodic boundary conditions. Namely the reference configuration is C, a circle of radius ^ (the 
points 0 and L are identified), 

or 

• Dirichlet boundary conditions. Namely u and its first n — 1 derivatives vanish at 0 and L. 

With above boundary conditions no boundary terms appear when performing the integrations by 
parts needed to obtain the equation of the motion (12.4p below. 

Note also that the maximal order of the spatial derivatives appearing in the equation of motion (12.4h 
is 2n. 

The equation of motion, as a consequence of the stationary action principle and the boundary condi¬ 
tions, is (with D^u = u) 

n 

u = - ^(-l)"T)“5^„$(u, Du,..., D^u). (2.4) 

Q( = 0 

We could also include, in the present context, a given external potential with a very minor effort. We 
avoid to do so for notational simplicity. 

Now we specify 4> by assuming that 

‘j’(i) = ^(e,Qe) + ^(i) ( 2 . 5 ) 

i.e. the quadratic part of $ is a quadratic form in terms of the displacement and its derivatives, contained 
in the vector C- Q = {Qa,p}a p=o ^ symmetric (without loss of generality) constant matrix with 
Qn,n A 0. 

On the non-linear part R we shall do suitable assumptions later on. We start by requiring that 

i?(0) = 0, R{Q = 0(\if), (2.6) 

namely the quadratic part of the interaction is fully expressed by the matrix Q. 

The fact that ‘h is not depending explicitly on x is consequence of the macroscopic homogeneity of 
the continuum (although it may be strongly inhomogeneous at microscopic scales). This implies that Q 
is constant. 
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As a first step we show that, in contrast with the fairly generality of the model, the quadratic part can 
be considerably simplified. Indeed, symmetrizing, integrating by parts and using the periodic or Dirichlet 
boundary conditions, we get: 


U = 


1 - ^ f 

2 


D^uDf^u 


a,/3=0 
2n 


-1)“+ (-!)<’] 


4 ^ '5'^- 

7=0 a,f}>0: 
a+l3='y 

1 ” f 

lE E + / iwk|" 


(2.7) 


7= o,/3>0: 
a+/3=27 


1 ” /■ 

o.—n '7/ 


where 


7=0 




( 2 . 8 ) 


a,^>0: 

a+/3=27 


Note that in the first step in (I2.7p we have used the symmetry of Qa,i3 and in the second step we used 
that (—1)“ + (—1)^ = 0 if a + /3 is odd. In the third step we have again integrated by parts. 

As a consequence of this analysis, without loss of generality, we can assume of the form 

1 "" 

+ (2.9) 

CK =0 

with An 7 ^ 0 and the equations of motion are 

n n 

u = - ^(-I)“A„A“u - Du,..., D^u), (2.10) 

q ;=0 a =0 

where A = denotes the Laplacian. Note that in the linear part only even derivatives are allowed. 

2.2 Formal discretization 

In view of the construction of the mechanical (microscopic) system with a finite number of degrees of 
freedom, we introduce a finite lattice of mesh e in I. The lattice points are {0, £,2e... , ke,... Ne} with 
the obvious condition Ne = L. When considering periodic boundary conditions we clearly identify 0 with 
eN. 

We associate to each lattice point a microscopic particle of unitary mass labelled by the index i G 
{0,..., A} and denote by Ui the displacement of the particle i from the reference position ie. The array 
Ue = discretized displacement field. 

The discretized Lagrangian takes the form 

N 


where 


Ce{Ue,Ue) = ^ ^ “ U{Ue), 

i=0 

(2.11) 

N . n 

E ^ [2 E 24«| (-DX)* 1' + R{{DeUe)^) , 

i=0 a=0 

(2.12) 
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a even 


where {DeUe) i = {{DeUe)i]a=Q, 


Here Df and , defined as 


= 




OL — 1 

Us, a odd. 


(2.13) 


{D+Us)i 


e 


{D Us)i 


^i—1 

e 


are the right and left discrete derivatives respectively and A^, defined by 


(2.14) 


(AsUs)i = (D+Us)i = {D^ Z1+ Us)i = ^{Ui+I + Ui-I - 2^*). (2.15) 

is the discrete Laplacian. 

To complete the above definitions we need to define the discrete derivatives at the boundary. For 
periodic boundary conditions it is enough to use the following convention: for any /c G Z, 


MiV+fc = Uk- (2.16) 

For Dirichlet boundary condition, we have to think of the first and last n particles frozen in their 
reference position. Hence we assume the constraints 

Ui = 0, i G {0,..., n — 1} U {A^ — n + 1,..., A'"}. (2-17) 

The equations of motion are 

dU 

Ui = Fi, Fi = (2.18) 

OUi 

with the index i running from 1 to A^ in the periodic case and on the set of i’s for which Ui is not 
constrained in the Dirichlet case. We notice that the choice of the right derivative (as well as any other 
possible discretization) is arbitrary. The only restriction that we have is the mechanical realizability (in 
principle) of this system. We are going to discuss this point in the next subsection. 

We finally remark that F) depends on Uj, with \i — j\ < n. However this is an almost local contribution 
because n is fixed and those Uj’s influencing m are at macroscopic distance 0{e). 


2.3 Realizable syntheses 

The aim of this subsection is to show that, at least in the simplest case of linear forces, the above 
introduced discrete system corresponds to a system of particles interacting via two-body forces of range 
not larger than n. Therefore, it can be realized by suitably assembling mechanical elements. 

Let us consider the linear system introduced in (j2.4p with R = 0 and its discrete counterpart (j2.18p . 
It can be checked that 

n 

Fi = - Y,{-lfAkA^,Us{xi). (2.19) 

k=0 

Therefore, the force acting on the particle i is expressed as a linear combination of discrete derivatives up 
to the order 2n. 

We want to show that F) can be interpreted as the result of the action of a system of linear pairwise 
forces with suitable range. More precisely, we want to And e-dependent coefficients kij such that 

F) = ^ kij{uj - Ui) (2.20) 

j 


6 








and hence 


U(u\,... ^ ^ ^ kij(ui Uj) . (2.21) 

*J=i 

We prove below that for any p, 

{A’’u). = Y;^Kfj{uj-u.), (2.22) 

3 

with Kf- other suitable constants. Once (12.221) is proved, we can conclude that (12.201) holds with 

n 

hj = (2.23) 

p=0 

Note that the constants ki^j are not necessarily all positive even if the are all positive. 

The constants Kfj are given by the recursive equation (|2.28p below. It implies that, for any p, Kf ■ 
vanishes for \i — j\ > p, thus kij = 0 if |)i — j| > n. Moreover, in the periodic case Kfj depends only on 
the difference i — j and is symmetric in the exchange i j and hence the action-reaction principle is 
satisfied. 

We prove (|2.22p by recurrence. 

For p = 1, we have 

— E T Ui—\ — E (u^-i-l U^) T E (u^—1 u^). (2.24) 

Thus p2.22p is verified with 

= Kli_i = and Kjj = 0 otherwise. (2.25) 

Suppose now that (|2.22p is true for p = £ — 1: 

(Af-iu)i = - Ui]. 

3 

Then, 


{Alu)t = = 5] Af-'KA.u), - (A,u).] 

3 

= E Ki7j^[^~‘^{'^3 + 1 - '^j) + - Uj) - E~‘^{Ui+i - Ui) - E~‘^{Ui-i - Ui)] 

3 

— ^ ^ ^ij (^J+1 '^i) ^ (^1 '^i) T ^ ^ i.'^j 

3 

- e“^(uj+i - Ui) - E~‘^{ui-i - Ui)]. (2.26) 

Using the change of index j + 1 —)• j in the first term and j — 1 ^ j in the second, we have 


i^iu)i = Y - Kj ^(“i “ + ^6+1^ ^(“i “ 


- e - Ui) - Ui) -E “^KYiui-i - Ui)]. (2.27) 

Thus, (12.221) is verified with the following recursive definition of Kfj: 


■ = E~‘^ 


K'-l. + K'ji - 2K'-‘ - (4+1., + 4_..j) 


(2.28) 


for ^ > 1 and given by p2.25l) . 

Equations p2(28p and (I2.23P solve definitely the posed problem of identifying the topology of the 
microstructure connections, since they provide the coefficients /cjj only in terms of the coefficients Ap 
that characterize the continuous formulation of the macroscopic description of the elastic problem. 
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3 A rigorous result of convergence 


In this section we prove a convergence result of the discrete model introduced in the previous section to 
the prescribed continuous systems in the limit as the scale parameter goes to 0. We show the convergence 
of the solution of the discrete system to the continuous one in the energy norm of the system. To clarify 
the argument without the use of cumbersome notation, we present first a paradigmatic case for which 
we discuss both periodic and Dirichlet boundary conditions. The more general case is considered in 
Subsection 3.2 where we give the convergence proof only in the periodic case although the argument can 
be straightforwardly extended to the Dirichlet boundary conditions as well. 

For the reader convenience we rewrite the Lagrangian we are going to consider in this Section, namely 

C{u^u) = ^ [ dx\u{x, t)\‘^ — - [ dx\D^u\‘^{x,t))'^ — f R{u, Du, D‘^u ...). (3.1) 

^Ji 2 Jj Ji 

As we shall see later on, we will con sider only nonlinear terms R depending on u and the first derivative 
only. 

3.1 The case - dynamic Euler-Bernoulli beam: “Elastica” 

3.1.1 Periodic boundary conditions 

We consider the Lagrangian (|3.1I) with Ag = = 0 and A 2 = 1. Moreover we focus on the linear case 

1? = 0. Thus we have the following linear initial value problem in the circle, C: 

d^u 2 

^ = (3-2) 

u{x,0) = uo{x), u{x,0) = vo{x). (3.3) 

It is well known that there exists a unique classical solution as the initial data are assumed sufficiently 
smooth. 

More precisely we assume that 

uq G , Vo G with s > 6,r > 4, (3.4) 

where denotes the Sobolev space endowed with norm 

S 

IIuIIhs = ^ ||D^u|||, 
e=o 


and II • lip is the U’{C)-noTm. 

In this way, by using the well known energy method, we can prove the propagation (in time) of the 

regularity for u and ii, yielding, in particular, u G C^{C) (as consequence of the obvious inequality 
Iklloo < C\\u\\hi). 

Next we consider the mechanical system of N particles, with coordinates Ui, i = 1,... N, whose 
Lagrangian is given by (I2.1ip again with Aq = Ai = 0, A 2 = 1 and R = 0. The equation of motion are 
explicitly 

Ui = ^(-Uj+2 + 4uj+i - 6ui - Ui-2 + dttj-i) i = l...N, (3.5) 

with the convention ujv+fc = Uk for any A: G Z. 

We want to compare the solutions of (|3.2I1 with the corresponding ones of (|3.5I) . To do this we first 
set 

Us{x,t) = Ui(t) if X G [fe, (f + I)e), fG {!,..., A}. (3.6) 
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In other words we introduce a function Ug which is the step, left continuous, function (constant in the 
lattice interval) taking the value of the nearest left point of the lattice. Problem (j3.5l) is rephrased 
accordingly: 

Ui;{x,t) = —A‘^Ui;{x,t) x^C, (3.7) 


where 


Asu{x) = Df u{x) 

= ±-(rt(x ± e) — u{x)). 


(3.8) 

(3.9) 


Notice that the Lagrangian (j2.11jl . with Aq = Ai = A 2 = 1 and R = 0, has the following continuous 
representation: 

(3.10) 




1 . . -.2 ^ / A / J .^^2 


C{Ue,U£) = / dx -ileixA) - T;{A£U£{x,t)y 


We suppose that, at the initial time, Ue,U£ are approximating u,u in the sense that 
Ue(x, 0) = no(ie), ^^(x, 0) = uo(ie) if x G [is, (i + l)e). 

Note that, by the conservation of the energy, we have 

1 


£[u{t)] := - / dx \u{t)\ +\Au{t)\ = £[u{0)], 

^ Jc 


as well as 


Se[u, 




dx 


\Us{t)\‘^ + \AeUe{t)\‘^ =^£^( 0 )] 


Next we introduce the following function which controls the deviation of from u: 

= ]- I dx {U£{x,t) - u{x,t)y + {Ue{x,t) - u{x,t)y + [A£(Ue(x,t) - u{x,t))f 

^ Jc 

Computing the time derivative and using the equation of motion we get 

wyt) = / dx{Ue{x,t) — u{x,t)){u£{x,t) — u{x,t) + U£{x,t) — u{x,t))+ 

Jc 

j dxA£(u£{x,t) - u{x,t))A£(u£{x,t) - u{x,t)) = 

+ / dx{u£{x,t) — u{x,t)){Ue{x,t) — u{x,t)) 

Jc 

— / dx{u£{x,t) — u{x,t))A^^{Ue{x,t) — u{x,t))\ 

Jc 

+ J dx{u£{x,t) — u{x,t)){A^u{x,t) — A^^u{x,t)) 

+ j dxA£{ue{x,t) - u{x,t))A£{u£{x,t) - u{x,t)). 

Now consider the following discrete integration by parts formula, namely 

[ f{x)Dfg{x) = - [ Dff{x)g{x) 

Jc Jc 

valid for any couple of bounded functions / and g. 


(3.11) 

(3.12) 

(3.13) 

(3.14) 


(3.15) 


(3.16) 
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If we apply the above formula twice we conclude that the second and fourth terms in (I, 1. 151) cancel 
each other. On the other hand the first term is bounded by 

- [ dx\u — + \u — Ue\‘^ < W. 

^ Jc 

The third term is bounded by 

^ J^dx{ueix,t) - u{x,t))‘^ + ^ Jj{A^ - Aj)u{x,t)\'^. 

Now the first term of the expression above is bounded by W. The second one, by the regularity of u and 
its derivatives up to the fifth order, is bounded, uniformly in x G C and in t in any bounded interval, 
by a constant cjg vanishing as e ^ 0. Here and in the rest of the paper G M denotes such a generic 
infinitesimal constant. 

In conclusion, by the Gronwall lemma, 


(3.17) 

so that We{t) is vanishing, because ITe(O) —0 by the regularity of u and the assumptions on initial data. 
We summarize above discussion in the following 

Theorem 3.1. Suppose that uq and vq satisfy (iTiD. Let u{t) he the solution to id. g|) and u^{t) be the 
step function defined by id.6\} with Ui{t), z = 1,... ,iV, solutions to lid.51) with initial data Ui{0) = uo{is) 
and Ui{0) = vo{ie). Then, for any t G M, 

lim Ws(t) = 0. 

£—^0 

3.1.2 Dirichlet boundary conditions 

For the Dirichlet boundary conditions we replace the circle C with the interval I = [0, L]. The equation 
(j3.2[) is well posed with the boundary conditions 

u{0,t) = u'{0,t) = u{L,t) = u'{L,t) = 0. (3.18) 

Remark: Several other boundary conditions may have an interest in engineering application and a 
physical meaning. For instance the conditions n(0) = n''(0) = 0, u{L) = u''{L) = 0, characterizes a 
beam with pivots applied at its endpoints, while the conditions which we considered here are relative to 
clamped-clamped beams. We do not consider in this paper the other possible boundary conditions, as 
the focus of this paper is different. 

Again, by using the energy method, we can construct solution with 77® regularity, by assuming 

uq G i7o n 77^, uo G 77q n 77^ with s > 6, r > 4. (3.19) 

Here 77 q (introduced to take into account the boundary conditions) is defined as the space of the 77^ 
functions vanishing in 0 and L, together with their first derivative. 

The corresponding discrete system is constituted by A"—3 particles with coordinates Ui, i = 2,..., N—2 

and 

uq = ui = un-1 = UN = 0 (3.20) 

are the constraints corresponding to the Dirichlet boundary conditions. 
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With this position, the explicit equations of motion are 

Ui = ^{-Ui+2 + - 6tti - Ui-2 + 4ui_i) i = 2...N -2, 

As before we introduce the left continuous step function 

Us{x, t) = Ui{t) if X G [ie, (i + l)e), i G {0,..., A" — 1}, 


(3.21) 


(3.22) 


but we find convenient to think of it as a function on M extended with value 0 outside I. Then (13.211) can 
be rewritten similarly to (13.7p as 


Ue(x, t) = —A^Us{x, t) X G le = (2e, L — e). 


(3.23) 


Note that the values of Ui are frozen for i = 0,1, A — 1, A, so that = 0 in We also think 

of the solution u of the continuous equation as extended with value 0 outside of I 
Next we introduce the function We{t) as 

^e{t) = 7 : I dx{ue{x,t) - u{x,t))‘^ + 1- I dx[Ae{ue{x,t)-u{x,t))f. (3.24) 

2 Jr 2 Jr 

Note that this function differs from the one defined by integrating on I instead of M because is non¬ 
local. It is actually larger and hence provides a stronger control of the convergence. Now we compute 
again the time derivative of W, as before and we get 

Wi;{t) = / dx{u^{x,t)—u{x,t)){ue{x,t)—u{x,t))+j dxA^{u^{x,t)—u{x,t))A^{u^{x,t)—u{x,t)). (3.25) 

By using twice the discrete integration by parts formula 


[ dxfDfg = - f gDf /, 

M t / M 


valid of any couple of bounded compactly supported functions / and g, the second term becomes, as 
before 

/ dx{us{x, t) — u{x, t))A‘l{us{x, t) — u{x, t)). 

Jr 

As for the the first term, we need to use the equations of motion ()3.2p for u and (13.231) for Wg. Note that 
the last ones hold only in Ig. Thus, using that tig = 0 in M — Ig, the first term becomes 

— / (ix('Ug(x,t) — m(x, t))(A^Mg(x,t) — A^ri(x,t)) — / ('Ug(x,t) — m(x, t))(—A^ m(x, t)) = 

Jle jR-/e 

— / dx{ue{x,t) — u{x,t)){A^^Ue{x,t) — A^u{x,t)) + / (-Ug(x,t) — t6(x,t))A^Ug(x,t)). (3.26) 

t/M i/M — /g: 

By adding and subtracting the term dx{ue{x, t)—u{x, t)){A^Us{x, t)—A^u{x, t)) the above term becomes 

— / (Jx('Ug(x,t) — n(x,t))(A^Ug(x,t) — Agtt(x,t)) — / dx{us{x,t) — u{x,t)){A'^u{x,t) — A'^u{x,t)) 

Jm. Jm 


+ 




(■Ug(x,t) - u{x,t))A‘^Ue{x,t)). (3.27) 


Putting together all these terms we conclude that 

W = — / dx{u£{x,t) — u{x,t)){A‘^u{x,t) — A‘^u{x,t)) + / (tig(x,t) — u(x, t))AgUg(x,t)) (3.28) 

IR.— le 
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The first term in the right hand side of (j3.28j) goes to 0 as in the periodic case, by the regularity u. The 
second term is the novelty of the Dirichlet case. In order to estimate it, note that hg = 0 outside of 
hence we need to estimate J^_j u{x, t)A‘^u^{x, t)). 

By the boundary conditions on u, {u = 0 and u' = 0 in 0 and L for any t), it results (by our assumptions 
\Au{x,t)\ is bounded) 

|'u(x, t)| <sup |Att(x, X G Ig. (3.29) 

2 xex 

Furthermore 

A^U£(x) = e~'^{AsUi;{x + e) + A^Usix - e) - 2AeUe{x)). 

Rewriting the total energy (13.1311 in a more explicit form, 

^ N-2 ^ N 

^[Ue] = 2 X] + 2'^A^eUe{ie)\^, (3.30) 

i=2 i=\ 

we obtain, at any time and for any x in X;., 


A£U£(x)| < 


V2Eq 


(3.31) 


where Eq = £{u{0)) is the energy of the initial data. Hence 

sup|AgUe(x)| < 4:\/2Eoe~^. (3.32) 

X&I 

Combining (j3.29l) and (|3.32ll and using the fact that the integration is restricted to the set I — I^, whose 
measure is 4e (remind that u = 0 outside X), we conclude that 


/ 


-le 


u{x,t)A‘^Ue{x,t)) <C\fe 


The rest of the argument proceeds as before and we conclude that Wi.{t) —)• 0. 
We summarize above discussion in the following 


Theorem 3.2. Suppose that uq and vq satisfy (I3.19h . Let u{t) be the classical solution to iS.SM with 
boundary conditions 113.18\) and initial values 113.3]) and Us(t) be the step function defined by /i3.22\) with 
Ui{t), i = 2,..., N — 2, solutions to i3.21]) with initial data Ui{0) = uofie) and 'Uj(O) = vo{ie). Then, for 
any t G M, 

lim We{t) = 0. 

£—^0 


3.2 A n-th gradient case 

Now we extend the previous argument to the more general setup corresponding to the Lagrangian (|3.1I) . 
restricting the discussion to the simpler case of periodic boundary conditions. The Dirichlet boundary 
conditions can be handled as in the previous subsection but we avoid here unnecessary complications. 
We assume the following conditions: 

1 . 

Aq > 0, An > 0, Aa >0, a = 1,..., n — 1 (3.33) 

2. We have already supposed that R(0) = 0 and R{f,) = 0(|^p). In addition we assume that, for 
n = 1, R depends only on u and, for n > 2, R depends only on u and Du. Moreover we assume 
R G C2’"+2(m2), 
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Remark 3.1: The positivity assumptions on the ^a’s with a = 1,... , n — 1, can be relaxed. In facts, let 
us define, for some eo > 0 , 


K = sup 


sup 


^ - '- j :- 11 p. 11 2 5 

eG(0,eo) Ji:||Dett|| 2 <l ll'^e'^112 

with the supremum on u taken on all u with 0 average. Then it is enough to assume 

n—l 


(3.34) 


E i-^-i 

a=l:Aa<0 




(3.35) 


to make the argument of the proof still working. This remark allows us to consider, for instance, the case 
il = (—— 7 A)u, with 7 sufficiently small, excluded by (|3.33p . 

Remark 3.2: The assumption on R concerning its dependence on u and Du only, is restrictive. We do not 
expect any surprise in assuming an explicit dependence on some higher derivatives. However, as we shall 
see in the course of the proof, more general assumptions would complicate the algebraic manipulations in 
dealing with the discrete derivatives in a consistent way. 

As regards the initial data we assume 

Uo G H 2 n+ 2 ^ ^ ^n+2^ (-3 

and, as before, the regularity is propagated. Clearly u G 


The explicit equation is 


u + ^(—l)"AaA"u + d^^R{u, Du) — D[d^j^R{u, Du)] = 0. 

a=0 

Note that thanks to the energy conservation, 

f 1 ” 

£[u]= / dx Aq,|T)'^uP| + R{u, Du), 

0=0 


(3.37) 


(3.38) 


we get immediately an a priori bound on the norm of u, ii and D'^w. 


u 


dx 


ftp + Aolttp + AAD'^u 




< ^[w(0)]. 


(3.39) 


Now we remind the discrete counterpart of the above setup, which corresponds to the discrete La- 
grangian (|2.1ip . Using the discontinuous function 


Us{x,t) = Ui{t) if X G [ie, (i + l)e), 
as in the previous section, the discrete Lagrangian can be written as 


Cr = dx 


L 


1 1 X \ 

-|Ue(x,t)P - - '^Aa\D^Ueix,t)\‘^ - R{Ueix,t),Dfue{x,t)) 
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a=0 


(3.40) 


(3.41) 


We can write the associated equations of motion in terms of u^ as 

n 

lie + ^(-1)"^oA"u£ + d(^QR{Ue,Dfu) - D~[d^^R{Ue, Df Ue)] = 0. 


(3.42) 


Q = 0 
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Also for the discrete system the energy conservation holds. Thus we have that 


f 1 ^ 

£e[Ue]= / dx + '^Aa\D'^u\‘^'^ + R{Ue,Dfue 

dC Q=o 


is conserved and hence, using that i? > 0, we have the inequality 


\L 


dx 


Us I + ^01 I + I De I 


|2 


< £e[Ue{^)]. 


(3.43) 


(3.44) 


Since > 0 and > 0, the existence, globally in time, for the solution to the discrete system follow 
from this bound. 

We start by proving the convergence of the discrere system to the continuous one in the linear case, 
namely when R = 0, 


u + y^(—l)"74Q,A"n = 0. 

a=0 

Similarly, the discrete system becomes 

n 

ii, + ^(-l)“^«AX = 0. 


(3.45) 


(3.46) 


a=0 


We introduce 
1 


We{t) = - dx\\u{x,t) - Ue{x,t)\‘^ + \u{x,t) -'Ue(a:,t)A + dxS^ Aa\D^[u{x,t) -Ue{x,t)]f. 

2 7c 7c ^ 

(3.47) 

The time derivative of W is: 


—We = / dx\ (ii - 


|(u - Ue)iu -Ue + U-Ue) + '^ AaD'^ {ii - Ue)D'^{u - U£)| 


ci=0 

n 


n 10 

/ dx{u - Ue)|(u -Ue + U-Ue) + y](-1)"A"(u - Ue)| 

0=0 

» n 

/ dx{u - Ue )\(u - Ue) + W(-1)“4 

7c 


Q :=0 
n 


A^Ue - A“n - Afue + A> 


} (3.48) 


n 10 

/ dx{u — Ue)^U — Ue + ( —1)"(A“u — A"u)| 


Q = 0 

In the second step we have integrated by parts a times, in the third we have used the equations of motion. 
The last step follows by canceling two equal terms with opposite sign. 

Hence the linear case goes exactly as in previous subsection, because | A“u — A"u| < uje for a < n and 
u e C2«+i(C). 

We summarize the results for the linear case in the following 

Theorem 3.3. Assume R = 0 and suppose that uq and vq satisfy (|3.36l) . Let u{t) be the classical solution 
to ^3.3T\) and Ue{t) be the step function defined by L3.6\} with ufit), i = 1,..., A, solutions to i3.45^ with 
initial data Uj(0) = uo{is) and Ui{0) = vo{is). 

Then, for any t E [0,T], 

\\u{t) — Ue{t)\\e —)• 0, as £ ^ 0, 


14 













where || • is the e-dependent norm defined by 




k=l 


(3.49) 


Next we consider the nonlinear case. Now the equations of motion are (I3.37P and ()3.42p for the 
continuous and discrete system respectively. Defining by (I3.47p . by the same computation, we have, 
again using the summation by parts formula, 




« n 

/ dx{u - Ue)^{u - + il - AaDf{u - Ue){Df{u - ufi) 

a=0 

» n 

/ dx{u - U£)|('U - Us + U - Ue) + - Ue)| 

C ^ —n 


- A^T/e + A^u 


+ 


^ u^e -r ' 


[ dx{u - Us){ {u - Us) + y~](-l)“A 

Jc ^ „=o 

— d^QR{u, Du) + Dd^-^^R{u, Du) + D^Us) — Dfid^^R{us, D^Us) | 

» n 

/ dx{u — U£)|u — Us + ^( —l)"^a(A"u — A"u) 

— d^QR{u, Du) + Dd^^R{u, Du) + d^^Rius, D^Us) — Dfid(^^R{us, D^Us) | 
To control the non-linear terms we proceed by estimating: 

Ti = d^^R{u, Du) - d^^R{us., D+ Us) = TI + Tf 


(3.50) 


+ 


and 

where 


T2 = Dfi[d^,R{us,Dtus)] - D[d^,R{u,Du)] = T^+D 
tI = d^QR{u, Du) — d^QR{u, D^u), 
rf = d^fjR{u, Dju) - d^^R{us, Df Us), 

Tl = D7[%i?(u,D+n)] - D[%i?(u,Dn)], 


and 


(3.51) 

(3.52) 

(3.53) 

(3.54) 

(3.55) 

(3.56) 


r| = Dfi[d^,R{us,Dtus)] - D,-[%i?(u,D+u)]. 

The bound ()3.39p and Poincare inequality imply that the L°° norms of tt, Du and D^u are bounded 
uniformly in £. Thus, by the local Lipschitz continuity of d^^R, we have 

|Tj^| < C\Du — D^u\ < ojs, 

by the regularity of u. Thus, by the energy bounds (I3.39P and (13.441) we get 

J dx{u-Us)Tl < C[£{u{0)) + £siUeiO))]^UJe, 

To control we need L°° bounds for Us and DsUs- They follow from the conservation of the energy 
for the discrete system by means of the following 
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Lemma 3.4. Let f be a step function on C left continuous in the points is. Suppose that 

\\f\\Hl= [ dx{\f\^ + \Dtff) 

® Jc 

is bounded. Then 

\\f\\oo<C\\f\\Hl. 

Proof. Let xq = ioe be any point such that |/(xo)P < Note that such a point does exist 

otherwise we would obtain a contradiction (f^ dx\f\‘^ > (ix|/p). For x = (zq + k)e we have 

k-l 

= f{xo) + J2[f{x(} + {h + l)e) - f{xo + he)]. 

h=0 


Since 


\f{x + e)-f{x)\ = \[f{x + e) + f{x)][f{x + e)-f{x)] = e[f{x + e) + f{x)]D^ f\ 

< + e) + fix)f + ie|D+/j^, 

we conclude that 

□ 


Lemma 13.41 and the energy bound (j3.44p imply that the L°° norms of and Df 'are bounded 
uniformly in e. Thus we can use the Lipschitz continuity of to get: 


L 


dx{u — Me)rf 


< KWe 


with K the Lipschitz constant of d(^^R in the ball of radius max{||u||oo, ||L>u||oo, ||'Ue||oo) H-DeUgHoo}- 

The bound of T 2 , involving discrete derivatives, requires the following chain rule formula for the 
discrete derivative of a composite function: 


Lemma 3.5. If f has continuous first derivative f, then for any function g and for any x there exist 
^e,x G (0,1) such that 

Dtfiaix)) = f'{Ce{x))Dfg{x), with Ce{x) = g{x) + eXe,xDfg{x) 


Proof. By the mean value theorem, for Df~ we have 

rgix+e) 

Dtfiaix)) = e~^[f{g{x + e) - f{gix)] = s~^ / dzf'(z) = e~^lg(x + e) - g(x)]f'(C) 

■tg(x) 

for a suitable ( in the interval with extremes g(x) and g(x + e): ( = g(x) + Xe,x[ 9 {x + e) — g{x)] = 
g{x) + £Xf;^xDf g{x) for some G (0,1). In the same way the statement for follows. □ 


By the chain rule. 


^^2 = dl^,e,idi{<fe{x),D'^u)D^ u - dl^^^^R{u,Du)Du + d‘^ 2 R{u,geix))AeU - df^R{u, Du)Au, 
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where 

Ceix) = u{x) + eXe,xD~u{x) 

and 

r]s{x) = D^u{x) + £^J.e,x^eU{x), 

with Xe^x G (0,1), ^Me,x G (0,1). But 

dlo,i,R{Ce{x),D+u)D-u - dl^^^R{u,Du)Du = 

Du)[D-u - Du] + D-u[dl^^^^R{Ceix), Df u) - Du)]. 

The smoothness of u and D^u and the Lipschitz continuity of ^^R yield 

\dl^,RiCeix),Dtu) - dl^^^Riu,Du)\ < C\Dtu - Du\) < u;,, 

so also this term goes to 0 by the regularity of u. 

Similarly, 

d‘^ 2 R{u, ri^{x))A^u — d^ 2 R{u, Du)Au 

= d‘^ 2 Riu, Du)[A^u — Au] + [d‘^ 2 Riu,r]e{x)) — d‘^ 2 Riu, Du)]A^u 

The first part goes to 0 by the regularity of u. By the boundedness of u, Du, Dfu and A^u, we can use 
the Lipschitz continuity of d%R to get the bound 

SI 

|5|2 -R(u(x),%(x),) - 5 |2 D(u(x),Du(x))| < K\rie{x) - Du{x)\. 

Since %(x) — Du{x) = D+u(x) — Du{x) + £^e,x^£u{x), 

\rie{x) - Du{x)\ < \Dju{x) - Du{x)\ + e^e,x\AeU{x)\, 

and hence 

|9|2i?(u(x),%(x)) - 9|2i?(u(x),Du(x))||A£u(x)| < (|D+u(x) - Du{x)\ + e/U£,a;|A£u(x)|)|A£u(x)| 

But 

|A£u(x)| < KAe - A)u(x)| + |Au(x)|. 

By the propagation of the initial regularity, ||Au(-,t)||oo is bounded for any t G (0,T). Therefore 
|9|2 -R(u(x),??£(x))A£u(x) - d'^ 2 R{u{x),Du{x))\\Ai;u{x)\ < uje. 

As for the term , we use again the chain rule: 


r| = D-[%i?(u„D+u,)] - D-[%i?(u,D+u)] = 

dl^^^^R{C,e{x),D'^Ue)D~Ue - dl^^^^R{Q{x), Dju)D~U + d'^2R{Ue,r]eix))AeUe - d‘^2R{u,fjeix))A£U, 


where 

(^(x) = Ueix) + £Xe^xD~Ue(x}, Ce(x) = u(x) + £X£^xD~u{x) 
%(x) = Djue{x) + £lle,xAeUe{x), fje{x) = D+u(x) + £He,xAeU{x). 
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We use the energy bound and Lemma [3l4l to get the boundedness of Q and Dfu^ and thus the Lipschitz 
continuity of (9|^ so that 

\dl^^^R{Q,,Dtue)D-Ue - dl^^^R{Q,Dtu)D-u\ 

< \dl^^^R{l,,Dtu)\\D-u, - D-u\+K\D-uMe{x) - C.(x)|. 

But 

lC£(a^) - C£(a:)l = \ue{x) - u(x)| + e{\D~Ue\ + \D~u\), 

so that 


L 


\dx\u - Ue||i:>+[%i?(ne,D+tte)] - [%D+u)]| < CW^ + '^e^{8[u{R)\ +£:£[Ue(0)]) 


The term d%R{u^^ri^{x))/S.i,Us — d^R^u, fj^{x))A^u is more delicate because, in order to use Lipschitz 

SI SI 

continuity, we need to bound ^(x) and hence the supremum of A^u^. But Lemma 13.41 and energy 
conservation are not enough when n = 2. However we need really to bound eA^u^ and we can take 
advantage of this extra e. Indeed, by the energy conservation and the positivity assumptions on R and 

eAn,n{D^,Ue{x)f < C, 

and hence 


^/e 


(3.57) 


implying that \rie^x\ ^ C (if n > 2 we get a better estimate). Thus we have 

\dfR{us,ris)AsUe - dfR{u,fjs)Aeu\ < \d‘^2Riu,fis{x))\\Ai;U - A^Usl 
+ \d‘^ 2 R{u,fjeix)) - d‘^ 2 R{u,r]^{x))\\AeUe\ + \d‘^ 2 R{ue,%ix))\ - d‘^ 2 R{u,r]eix))\\AeUe\ 

But, by (13.5711 . 

\r]six) - fis{x)\ = \D~{Ue{x) - u{x))\ + e{\A^Ue\ + lAeUl) < LVe- 
Therefore, by the Cauchy-Schwartz inequality and conservation of energy, we obtain that 

J dx\u - Us\\d‘^ 2 R{ue,rjs)AsUe - d‘^ 2 R{u,f]s)Asu\ <CIL + W£. 

Collecting all the terms, we conclude that there is are constant C > 0 such that 

d 


dt 


W,<CWe+U„ 


and hence, by Gronwall lemma 


\Ws{t)\ < |iy£(0)|e^‘ + fe^V as e ^ 0. 


We summarize the results in the following 

Theorem 3.6. Suppose that uq and vq satisfy (|3.36p . Let u{t) be the solution to |g.77| ) and u^it) be the 
step function defined by US.6\) with Ui{t), i = 1,..., N, solutions to \3.45^ with initial data Uj(0) = uo{ie) 
and Ui{0) = vo{ie). Then for any t > 0, 


\u{t) — u^{t)\\ —0, as e —)• 0, 


in the norm defined in \3.49). 
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